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Global,  n^wnil.  md  kxal  modds  of  die  Eaidt's  magnetic  field  have  a  wide  ran^  of  military 
anptkatkM  bodi  in  the  tradttiottal  warn  aa  aids  to  navigation  and  attitude/heading  reference 
lyaleaM  (AHRS)  for  aocraft,  imtailea.  tanka,  ahipa,  and  remtMdy  pUoted/operated  vehicles 
(RFV/ROV)  and  in  the  modem  sense  aa  tools  of  analysts  with  reflect  to  submarine  and 
HbicRiaean  deiection/conoeahiient,  dectronic  jamming,  endoatmospheric  dtar^  particle  beam 
wmponry,  and  An  moniloring  of  nuclear  test  shea.  Beyond  military  issues,  these  modds  also 
have  exteMtve  sciemilic  appKctfioBS.  Global  magnetic  fidd  modds  for  instance  are  used  as 
boundary  coaditioas  at  die  oore>aMUitk  interface  to  evduate  competing  dieories  of  die  Earth's 
magnelohythodynaraic  dynamo.  Regiond  and  local  fidd  modds  are  used  for  minerd  resoivce 
evahialion  and  ate  oriticd  to  the  undentandn^  of  the  gedogicd  history  of  die  Earth's  crust 

These  wv*d*t*  are  phenomenologicd  in  origin,  which  means  that  diey  require  nnvey  data  to 
support  them.  Furthermore,  smoe  the  Earth’s  Main  (Cote-generated)  magnetic  fkhl  varies 
slowly  but  cmiically  with  tane,  it  is  neccsswy  to  perform  globd  resurveyt  of  ^  fi^  at  teast 
pmiodkaQy.  if  not  oontnutoosly.  Low-altitade,  polar  orbiti^  satellite  platforms  provide  the 
moat  efBdeni  masm  adnevii^  die  required  ipi^  unifonnity  of  coverage  and  corresponding 
lempord  cohnsivenesa. 


Tha  Folar  Orbitng  Geosa^pietic  Survey  (POGS)  samllite  represents  the  Navy's  corttmuing 
allbrt.  aa  part  of  the  Prti^  MAGNET  program,  which  was  established  at  the  Navd 
Oceanoyaphic  Office  in  1951,  to  secure  the  necessary  data  to  maintain  accurate  magnetic  fidd 
10  support  aU  vahdaied  DOD  requirements  for  diese  models.  In  paiticalar,  this  data  is 
intended  to  support  dw  1995  Epoch  World  Mi^netic  Modd  (WMM-9S).  At  die  nme  time,  the 
vaiue  of  diis  data  for  scientific  purposes  is  folly  recognteed  Cdnsequendy,  at  widi  all  odier 
h^lh-altitiide  PlOfeci  MAGNET  airamft  data  imended  for  use  in  the  World  Magnetic  Modeling 
(WMM)  program.  POOS  data  is  rootindy  beo^  nade  available  to  the  general  public  through 
the  National  Geophyakal  Data  Center.  Boulder.  CO. 


T.E.  CALLAHAM 
Captain,  U.S.  Navy 
Commanding  Offker 
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The  spatial  and  taaporal  character  of  fully  calibrated  geoaagnetlc  data  collected  on 
aelacted  orbits  of  the  Polar  Orbiting  Ceoaegnetic  Survey  (FOGS)  satellite  •Isaion 
during  the  first  6  aontha  of  1991  is  analysed  with  respect  to  a  prelialnary,  degree 
12  apherical-haraonic  aagnetlc  field  aodel  derived  froa  aelacted  POGS  data  collected 
during  January  and  February  1991.  The  accuracy  of  global  aodala  derived  froa  these 
data  depends  on  the  vector  aagnetoacter  callbratioos,  themal  corrections,  satellite 
aagnetic  aapping,  and  the  tiae  syncing  of  the  satellite' a  quarts  clock  with  respect 
to  Che  aatallite’s  apheaeria.  which  is  referenced  to  Universal  TIsm  (UT)  and  which 
is  supplied  by  the  Defease  Napping  Agency's  tracking  network  (TRANET).  The 
calibration  and  correction  factors  applied  to  the  FOGS  data  set  ere  fully  docuaented. 
Additionally,  the  current  status  of  the  POGS  satellite  and  the  corresponding  data 
reduction  and  distribution  efforts  along  with  the  current  status  of  subsequent  POGS 
eaperiaents,  which  are  part  of  the  Defense  Meteorological  Satellite  ProgrM  (DKSP) 
Block  5  and  Block  6  satallita  aissions  known  as  DNSP/POGS  follow-on  expcriacots, 
are  presented. 


Magnetic  field  aodels,  local  field  aodels,  POGS  satellite. 
Project  MAGNET,  Defense  Meteorological  Satellite  Progrea  (DKSP) 
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1.  ffnmODUCTKM^ 


Tke  U.  S.  Navy*  Polv  Orbitiiig  Gcotmgckcac  Survey  (FOGS)  s^ilite  was  launched  1 1 
April  1990.  This  launch  cufamnaacd  the  efforts  of  aevoal  government  a^ncies  and  private 
contractors.  Prmci|de  uoa^  these  were  the  Naval  Oceanographic  OfRce  (NAVOCEANO);  the 
ConwMHdcr.  Naval  Oosanography  Comnand  (COMNAVOCEANCOM);  the  Defense  Mapf^ 
Agency  (DMA);  the  Air  Force  Space  Test  Program  (STP);  the  Navy  Space  Systcim  Activity 
(NSSA);  the  Offkx  of  Naval  Rmearch  (ONR);  the  National  Aeronautics  and  ^laoe 
AAntniitrabon  (NASA);  and  Detone  Systtms  lAcorporated  (DSl)  The  POGS  progrtun  is 
•poRSorod  by  COMNAVOCEANCOM,  which  faUs  under  the  auspices  of  the  Occanogr^ther  of 
the  hhivy  (N096).  It  »  NAVOCEANO^s  frsponsibittfy  to  operate  the  tatelliie  and  the  ground 
stations  which  codact  the  POGS  data,  while  n  ts  OMA‘t  respoMihility  to  trach  the  saielliie. 
asinf  hi  glohat  tradung  network  (TfLAHET),  asul  to  provide  the  taiciliie  ephemeris  which 
NA^^XTEANO  wergsa  whh  the  POGS  data.  The  integration  of  the  POGS  etperimem  with  the 
Adas  laanch  vshtcia  and  die  coordmaewn  with  the  OSl  contractors  that  built  die  satelltie  were 
the  leipotMSbtlify  of  STP,  widi  dm  aseMance  of  NSSA  and  (HfR  The  expenmeM  vector 
magnnifinnanr  was  bwili  by  NASA  and  was  cahbmed  md  tesmd  at  NASA's  Goddard  Space 
Cmm  (C^O.  which  opemaa  NASA's  Spacecraft  Msgnehc  Test  Faoltty  A  report  on 
dhia  fhciliiy  is  provided  by  HASAX^FC  (I9t4) 

The  POOS  teidhic  wae  Iwched  ntio  a  cncatar.  poinr  orbrt  as  part  of  a  multi-sateitite 
pra^act  known  as  P*t7  STACKSAT.  which  in  addMion  to  POGS  mcKidcd  two  commnmcatxMi 
lelePitBi  of  similar  etterior  desiyi  aa  indicaicd  m  F^arc  I  The  POGS  saieittic,  the  desi^ 
dbpacisrimica  of  whkh  are  haled  in  TaMe  1.  ts  at  die  top  of  the  stack  If  tos  52  soto  paneU. 
seme  of  which  mt  aied  aa  son  senson  to  provide  a  rough  esisnsm.  widim  a  tow  degrees,  of  the 
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T»blcl.  POGS  Dcs^  Oiuractnlstict 


•  EXPCIUMENT  MAGNETOMETER 
<  MOMFIED  DMSP/SSM  TRIAXIAL  FLUXGATE  (NASA) 

•  RANGE:  ±iSS3$«T 

-  RESOLUTION:  ±2aT 

-  LOCATH)N:  AT  TIP  OP  f-FT  CWB«  (NON<MAGNETlC)  BOOM 


•  ArmUDE  CONTROL 

•  GRAVITY  GRADIENT  STABILIZATION 

•  2B>PTBOOM 

•  1I1>LB  TIP  MASS 

•  4  HYSTERESIS  RODS 

-  '*Z’' TORQITNG  COIL 

•  45  AMP’TURN/METER 

«  PERPORMS  SATELLITE  PUP  MANEUVERS 


•  ATTITUDE  DETERMINATION 

-  ONBOARD  TRIAXIAL  PLUXGATE  MAGNETOMETER 

•  RESOLUTION:  ±5U«T 

•  LOCATION:  ONBOARD  SPACECRAFT 

•  f  SUN  nNSORS  LOCATED  AROUND  SATELLITE 

•  RESOLUTION:  ±  I* 

•  DETBIDICINBS:  AZ2MUTN  B  ELEVATION  (RELATIVE  TO  SUN-UNE) 


•  POSmONlNC 

•  RESOLUTION:  75  METERS  SPHERICAL  PROBABLE  ERROR  (SEP) 

•  GRAVITY  MODEL:  CEM-ttB 
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satellite's  attitude  (i.e.,  azimuth  and  elevation  relative  to  the  Sun-line).  The  satellite  is 
gravity-gradient  stabilized  and  is  tracked  by  a  global  network  which  monitors  the  satellite's  two 
navigation  beacons,  which  transmit  signals  at  separate  frequencies  at  20-second  intervals.  The 
satellite  carries  two  vector  magnetometers,  the  coarse  onboard  Navigation  magnetometer,  which 
is  used  for  attitude  sensing  of  the  satellite,  and  the  high-resolution  Experiment  magnetometer, 
which  is  located  at  the  tip  of  an  8-ft  nonmagnetic  boom.  These  are  oriented  relative  to  the  earth 
and  with  respect  to  each  other  as  indicated  in  Figure  2. 

Four  significant  problems  were  encountered  during  the  satellite's  initial  operations,  all  of 
which  have  been  overcome.  First,  the  satellite  was  launched  into  orbit  in  an  inverted  mode. 
There  was  a  50  percent  probability  that  this  would  occur.  As  a  consequence,  the  satellite's 
communications  antennae  pointed  spaceward  rather  than  earthward.  This  caused 
communication  problems  between  the  satellite  and  the  ground  stations.  Although  torquing  coils 
were  available  to  flip  the  satellite,  this  was  considered  a  risky  procedure  which  could  have 
snapped  the  gravity-gradient  boom  or  the  Experiment  magnetometer  boom.  Consequently, 
several  months  were  spent  redesigning  the  ground  station  antennae  patterns.  The  Master  ground 
station  is  located  at  Stennis  Space  Center,  MS.  A  second  ground  station  was  operated  by  the 
U.S.  Geological  Survey  through  1992.  It  was  located  at  Fairbanks,  AK.  This  ground  station  has 
since  been  replaced  with  another  ground  station  operated  by  the  U.S.  Navy  at  Keflavik,  Iceland. 
The  communications  problem  was  resolved  by  November  1990.  A  second  problem  concerned 
the  satellite’s  onboard  data  packing  and  storage  program.  The  Experiment  magnetometer's 
data  collection  syston  was  designed  to  store  1 -second  data  samples  in  10-second  packets.  The 
ftnrt  recmd  in  each  packet  was  a  complete  record,  while  the  following  nine  records  were 
l-secmd  ddta  valuM.  Unfortunately,  the  program  placed  the  Y-component  deltas  into  both  the 
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Y-«oa|iOQent  and  Z-comfionent  storage  slots.  This  crnN’  was  uncorrectable,  resulting  in  an 
effective  <kta  scnaple  rate  of  10  saaq;)les  po*  second  instead  of  1  sample  per  second.  The  thiid 
proMem  cooMmed  the  inal^ty  to  montUv  the  drift  of  the  satellite's  onboard  quartz  clock  vu 
adiediiled  timing  polsn.  This  problon  was  resolved  by  another  modification  to  die  ground 
atatkm  hardware  in  Oecendier  1990.  As  a  residt,  useful  POGS  data  begins  in  January  of  1991. 
Satdhte  operations  were  contimioas  from  that  time  (except  for  some  minor  down  time  due  to 
radimion  htta)  until  September  1992,  when  radiation  caused  the  fourth  {uoblem,  a  disruimcMi  of 
that  part  of  the  satellite’s  memory  which  contained  data  downlink  instiuctitms.  Omununicatitm 
with  die  satdlite  was  lost  for  5  months  but  was  reestablished  in  January  1993.  Frmn  diat  time 
forward  the  satellite  has  continued  to  operate  normally.  The  satellite's  orbit  characteristics  and 
other  minton  details  are  listed  in  Table  2. 

2.  MAGNETOMETER  CAUBRATIONS 

The  Experiment  magnetometer  was  calibrated  at  NASA/GSFCs  Magnetic  Coil  Facility, 
where  die  magnetometer's  scale  fiactors,  biases,  orthogonality  fturtort,  and  temperature 
characteristics  were  determined.  Scale  factors  for  the  Navigation  magneunneter  were  also 
determined  this  fsdlity.  Some  time  later,  when  the  Experiment  magnetometer  was  attached 
to  the  satellite  and  the  boom  extended,  a  static  magnetic-field  measurement  was  performed  to 
detenuine  die  satellite's  magnetic  signature  at  the  site  of  the  Experiment  magnetometer.  This 
test  was  performed  with  all  saseUite  electronic  packages  turned  on. 
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Table  2.  POGS  Minloa  Detafb 


•  ORBITAL  PARAMETERS 

•  APPROXIMATE  PERIOD: 

•  SEMI-MAJOR  AXIS: 

-  SEMI-MINOR  AXIS: 

-  APOGEE  RADIUS: 

-  PERIGEE  RADIUS: 

•  INCLINATION: 

•  MEAN  ALTITUDE: 


•  LAUNCH  PARAMETERS 

-  LAUNCH  DATE: 

-  LAUNCH  TIME: 

•  LAUNCH  SITE: 

-  LAUNCH  AZIMUTH: 


W.MMia 
7059.44  KM 
7059 KM 
7115.95  KM 
7002.94  KM 
09.8SDct 
4nJ3KM 


11  APRIL  1990 
0200  UT 

VANDENBERG  AFB,  CALIFORNIA 
100  Deg 


•  MISSION  OPERATIONS 

-  MISSION  LIFE  EXPECTANCY:  1  T03  YEARS 

-  TRACKING  ORGANIZATION:  DEFENSE  MAPPING  AGENCY 

-  TRACKING  BEACON  FREQUENCIES:  150/400  MHz 

-  DOWNLINK  GROUND  STATIONS:  NAVOCEANO,  SSQ  MS 

USGS,  FAIRBANKS,  AK 

U.S.  NAVY,  KEFLAVDC,  ICELAND 
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2.1  Tke  Navtgatkm  Magni^neter 

The  Navigali<m  vector  magnetometer  is  intended  to  measure  the  ambient  magnetic  field  to  a 
reacdution  of  only  ±511  nanoTnlas  (nT).  The  intended  purpose  of  this  magnetometer  is  to 
detennine  the  satellite's  orientation  primarily  just  after  orbit  insertion  and  satellite  spin-down. 
Eadi  vector  axis  has  its  own  scaling  fonnula  which  converts  from  counts  (i.e.,  relative  voltage) 
to  nanoTealas.  The  dimensi<xileaa  scale  factors  for  the  X-.  Y-,  and  Z-axes  of  this 
magnetometer  are: 


S.  «  0.92761 

(la) 

Sy  «  1.05508 

(lb) 

S,  -  0.99245 

(Ic) 

With  n,  ,  iiy  ,  and  n,  being  the  measured  counts  (range:  0  to  255)  from  the  Navigation 

magnetometer,  the  magnetic  field  for  this  magnetometer  is  reconstructed  for  each  vector  axis  in 

nanoTealas  as: 

B,  «  (n,  >  132)(60000.0/127)S. 

(2a) 

By  «  (ny  -  129)(60000.0/127)Sy 

(2b) 

B.  «  -(D,  -  130)(60000.0/127)S, 

(2c) 

At  dM  tiine  that  measurements  ate  taken,  the  axes  of  the  Navigation  magnetometer  (i.e.,  the 
onbotfd  magnetmneter)  are  orioited  as  indicated  in  Figure  2,  which  is  <  tt-handed  syston. 
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The  Navigation  magnetoaneto^s  X-  and  Y<axea  are  parallel  to  the  q>acecraft  coordinate  axes, 
ii^ii]e  the  Z<axis  is  anti>parallel  to  the  corresponding  Z-axis  of  the  jq>acecraft,  vtMch  is  the  ^>in 
axis.  The  leading  negative  sign  in  eq.  (2c)  inv^  the  Z-axis  of  the  Navigation  magnetometer, 
thus  making  it  parallel  to  the  q>acecraft  coordinate  Z-axis  and  in  the  process  convi?rts  the 
left-handed  Navigatitm  magnetometer  coordinate  system  into  a  right-handed  coordinai  lem, 
with  the  Z-axis  pointing  toward  the  earth. 

2.2  The  Experiment  Magnctoaicter 

There  are  four  aspects  to  the  calibration  of  the  high-resolution  Experiment  magnetometer. 
The  first  aspect  is  concerned  with  determining  the  calibration  factors,  \riiich  convert  die  sensor 
measurements  from  counts  (essentially  voltage  measurements)  to  nanoTeslas.  The  second  aspect 
is  concerned  with  the  precise  determination  of  the  orthogonality  (i.e.,  alignment)  of  the  vector 
axes  of  the  Experiment  magneUmieter.  The  third  factor  involves  the  determination  of  die 
magnetic  noise  generated  by  the  qiacecraft  and  sensed  by  the  Experiment  magnetometer  at  die 
tip  of  the  fuUy  extended  nu^nebmieter  boom,  when  all  electitmic  packages  are  operating.  The 
fourth  facUv  is  to  determine  the  toiqperature  senativity  of  the  three  axes  of  the  Experiment 
magnriometer. 

The  Experiment  magnetometer  is  of  the  same  type  that  was  designed  and  built  by  NASA  for 
die  Defense  Meteorological  Satellite  Program  (DMSP).  The  only  modification  in  the  design  was 
to  add  extra  dwrmal  mass  in  order  to  reduce  the  magnetometer’s  thermal  drift  to  something  less 
than  SO  nT/yr.  The  precise  thermal  drift  is  unknown.  However,  it  is  expected  that  this  drift 
should  be  die  greatest  just  after  the  satellite's  initial  insertion  into  orbit  and  that  after  a  few  days 
or  wedcs  would  setde  down  to  a  rather  small  value.  There  is  no  absolute  scalar  magnetometer  to 
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monhor  thu  drift  rate  and  no  odwr  direct  means  to  determine  it.  The  initial  analysis  of  the  dau 
after  neaidy  a  year  in  orbit  indicates  that  the  drift  is  much  less  dun  SO  nT/yr,  when  con^pared  to 
predictioos  of  1990  ^wdt  ^dterical-harmcmic  magnetic  fidd  modds  that  were  generated  from 
<teta  independent  of  POGS. 

The  calibratian  constants  for  die  POGS  Experiment  magnetometer  are  given  in  Table  3.  The 
"measured*  nugnetic  ftdd  vector  coo^ionents  based  on  these  calibradm  omstants  are 
computed  according  to  the  fdlowing  fonmilae: 


s 

B»  •  a*o  +  aiabia(«MMi)  Kt(b((aB4  ■“  bzo)  (5a) 

■a| 


s 

By  •  ayO  +  2|  ay  by»(iii— )  +  Ky  (by0bi^  "*  byo)  (5b) 

•-I 


B. 


a,o  £  a..b  ^  KaCbaftaa)  bao) 


(5c) 


The  ftetors  b.^^.  -  j  ,  b^,^,  ,  arui  b^,^,  ate  the  COARSE  bias  step  intervals  of  die 
measured  magnetic  fidd,  sriiile  factors  b^, .  b,^, ,  and  b^,  are  the  FINE  adjustments  of  the 
measured  magnetic  fidd.  The  division  of  the  magnetic  fidd  measurements  into  COARSE  and 
FINE  segments  is  rdated  to  die  compact,  binary-packing  scheme  used  to  omqxess  and  stme  die 
data.  Notice  again  diat  in  die  reconstruction  of  the  magnetic  fidd  vector  exm^pments  a 
coordinate  transformation  must  also  been  performed.  As  indicated  in  Figure  2,  du  X-axis  of  the 
Experiment  magnetometer  (Experiment  sensor)  is  aligned  along  du  ^lin  axis  of  die  satdlite, 
die  Z-axis  of  dus  magnetometer  is  in  the  X-Y  plane  of  the  spaceoaft's  coordinate  system. 
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TiUric  3.  POGS  Expcrinicat  Mj^^omctcr  Calibratfoa  ConsUatt 


^^^^.^^MPONENT 

CONSTANT'\,^^ 

X’AXIS 

Y-AXIS 

Z-AXIS 

K  |aT/cMatl 

1.99S240 

1.999900 

1.999344 

••  I■T1 

•M3SSJ6 

•66317.94 

-66379.24 

Ian 

+66319J7 

+66390.95 

a,  laTI 

+33196.44 

+33159.64 

+33199.73 

a,  |aT| 

+16599.49 

+16590.97 

+16593.73 

*4 

+9299.79 

+9291.02 

+9299.54 

«,  (■Tl 

+4149.44 

+4145.13 

+4147.95 

b,  {coaatal 

+2052.90 

+2056.90 

+2050.20 

n 


E<ptatu»<  (3«)  and  (3c)  tnuufonn  the  Experiment  magnetmneter’s  comdinate  axes  so  that  its 
Z-axia  is  pvaUd  to  die  spacecraft’s  spin  axis,  ^tdiile  die  magnetometer’s  X-axis  lies  in  the 
spacecraft’s  X-Y  plane.  After  this  transformation,  the  Experiment  magn^otneta  axes  are 
spproximatdy  aligned  to  the  coordinate  axes  of  the  spacecraft.  They  are  not  exaedy  aligned 
because  die  boom  upon  which  this  magnetometer  sits  is  a  flexible  helical  coil  which  can  twist 
and  bend  <hie  to  thennal  action,  satellite  motion,  and  initial  dqiloyment  orientation.  None  of 
diese  factors  are  monitored.  Consequoidy,  the  precise  orientation  of  the  Experiment 
magnetometer  rdadve  to  die  spacecraft’s  coordinates  is  unknown.  Coarse  estimates  of  this 
magnetometer's  orientation  can  be  made  by  coa^wring  its  vector  magnetic  measurements  to 
those  of  die  less  accurate  Navigation  magnetometer,  whidi  is  fixed  in  space  relative  to  dM 
spacecraft's  ooordinines. 

The  magnetic  field  values  derived  in  eqs.  (3a),  (3b),  and  (3c)  are  referred  to  as  die 
"Measured”  vector  components  of  the  earth's  magnetic  field.  However,  the  vector  axes  of  die 
Experiment  magnctomeiCT  are  not  precisely  Signed  perpendicular  to  each  other.  The  "True" 
magnetic  fidd  vector  components  associated  with  an  orthogonal  set  of  axes  are  obtained  from 
die  "Measured"  magnetic  field  components  via  the  following  alignment  «»Tection 
transformations 

(4) 

where  the  3  x  3  mmsformation  matrix,  the  elements  of  which  are  dimensionless,  is: 
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l.OOOO  -7.ti96xlO'*  -3.2410x10-'+  J  Yi.B; 

»i 

5 

-3.OI90xlO-»  0  9999  1  2820x  lO"*  +  J  T*.B; 


•  6150x10-*  +  2  9630x10- 


1.0000 


a«l 


(5) 


Tlw  croM  comtaott  we  listed  in  Tal^  4.  Fuither  details  coooeniing  these 

IH4SP>type  magnetoineter  calibrationt,  with  respect  to  the  determination  of  magnetometer 
calilintiott  coeffidentfs.  scale  factors,  and  ordiogooality  factors,  are  discussed  by  Chomay 
(IW7). 

The  third  stq>  in  the  caUfaratioa  procedure  is  to  correct  for  ambient  temperature  fluctuatkms. 
POOS  has  two  semperature  sensors,  one  to  monitor  the  msgnetoaieter*s  dectronics  package, 
dencned  as  T| ,  and  one  to  monitor  the  magnetic  aenaor,  denoted  as  T, .  The  former  is  onboard 
dw  spacecraft  and  latter  is  at  the  tip  of  the  boom  and  therefore  will  always  measure  slightly 
oohkr  temperatures.  The  temperanve  correction  factors  are  referenced  to  room  temperature 
(2(f  C).  Eadt  axis  of  die  magnrtomrtcr  has  its  own  correction  factor.  The  temperature 
corrections  applied  to  using  the  sensor  temperature  T,  are: 


B,  *  [1.0  +  0.0000«0(T>  -  20)]B.(n^  (6a) 

B,  «  [1.0  +  0.000046(Tt  -  20)]B,(ita4  (6b) 

B.  »  [1.0  +  0.000047(Ts  -  20)]B.(iv^  +  7.0  (6c) 
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TaMt4w  POGS  Eipcrtecat  MagMtMictcr  CrMS  Coupling  CoBstanU 


•  Tta  Ti.  *6. 


1 

2 

3 

4 


-3,474 1 19* 
•4.5S4 1 10*“ 

+2^24 1  ir“ 

+3.070  K  ir" 


-1.122  t  ir* 
-3.040 1  l(r“ 
+7.417  1 10*“ 
-1 J03  » 10*” 


+4.527 1 10^ 
-7.044  X  l(r‘* 
-3.747  X  10^* 
+4.050  X  10*" 


5 


•4.433  xir" 


-1 J17  X  10^ 


+7J77  X  10*“ 


The  UMt  conectioii  to  the  Experiinent  ma^teUMneto^  obswaticms  accounts  fcNr  the  magiMtic 
fidd  generated  hy  die  body  of  the  spacecraft  and  aoued  at  the  tip  of  the  magneumieter  botm 
Laboratory  measurements  found  that  the  only  correction  required  was  the  additim  of  a  7 
nanoTeala  bias  to  the  Z-t^n^ionait  of  die  tcnqimture-corrected  field.  This  bias  has  been 
included  as  put  of  the  Z-coeoponent  tempoature  correction  of  eq.  (6c).  This  bias  does  not 
inchide  corrections  for  the  radiative  effects  of  die  two  navigation  tracking  beacons  whidi  are 
turned  on  for  aboin  20  seconds  of  each  minute.  The  magnetic  field  sensed  at  die  tip  of  the  boom 
due  to  these  beacons  is  on  the  order  of  3  or  4  nanoTeslas.  However,  no  corrections  for  die 
magnetic  effect  of  the  beacon  transmissions  were  applied  to  the  data.  Instead,  a  flag  is  set  when 
the  beacons  are  on.  Then,  for  modeling  purposes  only  unflagged  data  ate  selected. 

3.d  SUN  SENSOR  ATTITUDE 

The  POOS  satelliie  has  S2  solar  panels.  Some  of  these  panels  have  been  configured  as  sun 
sensors  in  order  to  roughly  determtiie  the  satellitea  attitude  relative  to  theSun'line.  The 
parameters  measured  are  die  Azimuth  and  Elevation  of  die  spacecraft  relative  to  the  Sun-line. 
The  accuracy  of  these  pwameters  is  not  known,  but  they  are  certainly  no  more  accuntt  than  ±  3* 
under  the  beat  of  drcumsianccs.  The  accuracy  degrades  as  the  earth's  dawn-dusk  meridian  is 
approached  from  die  sun-side.  The  +  Z  axis  of  the  sun  sensor  coordinate  system  is  ptfallel  to 
die  ^laceoraft's  spin  axis,  which,  because  the  satellite  wu  inserted  into  Mbit  upside-down,  is 
Eerth-Pmntii^  However,  the  X-axis  and  the  Y-axis  of  the  sun  sensor  coordinate  system  are 
rotated  1 1.25*  counterclockwise  with  respect  to  the  corresponding  ^Uteccraft  axes  as  indicated 
in  Figure  3.  In  this  figure,  die  satellite's  axis,  which  is  Earth-Pointing,  is  directed  out  of  the 
|M^  Simply  used  as  nin  sensors,  the  solar  parul  data  are  valuable  ui  <fetenntning  when  the 
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POGS/SS/t  SA  TELUTE 


latcrior  SdMOMtk  •!  POGS  Satdiltc:  S«a  Scmior  X-Y  Axes  are  Routed 
IKZS*  Cwterciocfcwfat  wMi  Respect  to  Spacecraft  X-Y  Axes 


m^te  is  on  the  sun-ade  of  the  esrth  snd  what  it  is  in  t}»  earth's  shadow.  This  pmnits  easy 
iftterpeetalion  of  the  magnetic  field  re^wnses  al<mg  each  orbit  as  being  of  cru^  origin  or  of 
extcnial  kmosphearic  or  magnetoq^ieric  origin. 


4.  SATELLITE  POSITIONING 

There  nre  two  tracking  bncons  on  the  POGS  satellite.  Each  beacon  has  its  own  transmitter 
freqooicy.  The  beacons  are  tracked  by  the  TRANET  global  network  of  tracking  stations,  the 
d^  frosn  which  are  siqi^lied  to  DMA,  which  in  turn  processes  the  tracking  data  using 
sophinicrted  orbit  determination  programs  developed  by  the  Naval  Surface  Weapons  Center's 
Dahlgroi  Labonttory  in  Oahlgren,  VA.  The  POGS  satellite  trackmg  data  processed  by  DMA 
«e  refoenced  to  the  GEM«10B  gravity  model.  Numerous  tidal,  ionoqiheric,  and  atmospheric 
corrections  are  applied  to  die  tracking  data.  The  final  ephemeris  (i.e.,  orbit  determination)  is 
supplied  to  NAVOCEAhK)  at  1*010^16  tiroe  intervals  in  EardHCcntered,  Eafth>Fixed 
coorduiam.  relative  to  the  coorduune  system  of  the  GEM>IOB  ellipst^  with  sqiherical 
probable  oror  (SEP)  of  less  dian  75  mtfers.  Further  details  oonceming  satellite  oibit 
detennhwiion  are  discussed  by  James  W.  OToole  (1976)  and  by  CappeHari  el  al,  (1976). 
Relevant  mformation  omccniti^  geodetic  positioniim  is  given  by  Mudler  and  Eichhorn  (1977), 
while  relevant  information  concemom  earth  rotation  is  given  by  Moritz  and  Mueller  ( 1987). 

Due  to  die  recent  completion  of  the  Global  Posttioning  System  (GPS)  network  of  saidlites, 
DMA  hm  made  the  decision  to  disaiMide  its  ground-based  global  tteellite-tracking  network  and 
its  data  procesatng  support  group  bcginnnm  in  October  1993.  Consequendy,  die  status  of  the 
POOS  MteQHe,  which  is  sttll  espected  to  be  opentionsl  st  that  time,  then  conies  into  question. 
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There  are  altentadve  sources  for  this  tracking  data.  The  possibility  of  acquiring  ephemeris 
data  from  these  alternative  sources  is  now  being  investigated. 


5.  POGS  CLOCK  CALIBRATIONS  AND  DATA  TIME  SYNCEVG 

The  quartz-crystal  clock  on  the  POGS  satellite  is  accurate  to  about  one  part  in  10^  and  has  a 
Ribstantial  linear  drift  which  is  occasionally  interrupted  by  discontinuous  jumps.  The  drift  and 
discontinuous  JunqM  in  the  POGS  clock  readout  can  accumulate  to  a  timing  error  of  several 
seconds  over  the  course  of  several  months.  Consequently,  it  is  necessary  to  periodically 
caUtmue  the  POGS  clock. 

The  POGS  quartz-clock  calitnadon  procedure  involves  the  transmission  of  a  20-second 
timing-piilae  from  the  spacecraft  to  the  Master  ground  station  located  at  NAVOCEANO,  Stennis 
Space  Center.  MS.  The  start-time  and  end-time  of  each  pulse  are  known  a  priori  from 
pre-determined  schedules  telemetered  to  the  satellite.  Due  to  the  fact  that  the  onboard  computer 
may  be  performing  other  tasks  when  a  particular  timing-pulse  is  intended  to  begin,  the 
tmung-pulse  cuy  not  begin  exactly  when  scheduled.  In  this  event,  the  timing-pulse  is 
KitosnaticaUy  shortened  so  that  it  ends  exactly  20  seconds  after  its  intended  start-time.  Thus, 
measuriiqt  the  time  that  the  trailing  edge  of  the  timing-pulse  reaches  the  Master  ground  station 
using  a  taMium  time  standard,  making  cmrections  for  the  pulse  travel-time  from  the  satellite  to 
die  grtmnd  stirtion,  and  subtracting  the  intended  end-time  of  the  timing-pulse,  yield  the  time 
dtlfcratce  between  the  POGS  clock  and  Naval  Observatory  (Universal)  time  to  which  the 
ceaiiim  standard  is  adibrated.  The  Master  grouiul  station’s  cenum  standard  clock  is  accurate  to 
within  a  few  nuaroaeconds  per  year  and  is  dwrefore  quite  stable. 
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The  travd-tinie  correction  computation  requires  knowledge  of  the  slant-range  between 
die  satdlite  and  the  ground  station.  The  slant-range  is  detennined  from  an  orbit  computation 
algoarithm  developed  by  DSI,  the  prime  contractor  responsible  for  building  the  POGS  satdlite 
and  for  writing  the  satellite  and  ground  statimi  software.  The  orbit  parameters  for  this  algorithm 
are  regulariy  updated  using  values  supplied  by  the  Air  Force's  Satellite  Tracking  Facility  at 
Cheyenne  Mountain,  CO.  With  the  slant-range  given  in  nautical  miles,  the  travel-time 
correction  in  seconds  is  given  by  the  following  formula: 

T(tiavd)  =  6.183  X  10"*  R(aMi)  (7) 

The  time  difTerences  between  the  POGS  clock  and  the  cesium  standard  clock  for  the  period 
fttun  January  1991  to  June  1993  are  plotted  in  Figure  4.  Their  actual  numerical  values  are  listed 
in  Table  5.  Figure  4  indicates  that  the  POGS  clock  generally  runs  a  few  seconds  ahead  of  die 
cesium  standard  clock.  Disjoint  linear  drift  segments  ranging  in  duration  from  a  few  days  to 
more  than  a  month  in  length  are  clearly  evident  Occasionally,  the  computed  time  differences 
are  significandy  in  error.  These  errors  can  be  attributed  to  unusual  iono^heric  or  atmo^heric 
coi^idons  whidi  occurred  during  the  timing  pulse  and  also  to  human  error.  These  mors 
produce  the  occasional  but  obvious  anomalies  seen  in  Figure  4  and  are  discarded  beftne 
determining  the  POGS  clock  corrections. 

The  POGS  time  correction  procedure  involves  the  identification  of  the  beginning  and  end 
times  eadi  time  sub-interval  for  which  the  clock  drift  is  linear  and  the  determination  of  the 
clock-drift  rate  for  eiKh  of  these  sub-intervals.  The  time  tags  attached  to  m^netics  data 
acquired  during  a  particular  time  nib-interval  are  then  appropriately  adjusted  to  coincide  with 
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CALENDAR  DAYS  FROM  JANUARY 


Table  5.  Ttining  Pulse  DifTerences  (POGS  -  UT);  Units:  Seconds 


YEAR 

DAY 

DIFFERENCE 

YEAR 

DAY 

DIFFERENCE 

1991 

015 

+2.533 

1991 

241 

+9.234 

1991 

030 

+2.295 

1991 

246 

+9.549 

1991 

032 

+1.833 

1991 

256 

+10.153 

1991 

037 

+1.760 

1991 

259 

+10.283 

1991 

043 

-0J19 

1991 

263 

+10.549 

1991 

052 

+0.776 

1991 

266 

+10.711 

1991 

079 

-0.183 

1991 

273 

+10.981 

1991 

129 

+1.839 

1991 

277 

+10.736 

1991 

140 

+2.497 

1991 

281 

+10.925 

1991 

150 

+3J82 

1991 

284 

+11.070 

1991 

169 

+4.178 

1991 

294 

+9.852 

1991 

171 

+4J24 

1991 

297 

+9.794 

1991 

175 

+4.621 

1991 

304 

+8.642 

1991 

196 

+5.139 

1991 

308 

+8.830 

1991 

200 

+5.426 

1991 

308 

+8.830 

1991 

211 

+6.218 

1991 

322 

+9.571 

1991 

225 

+8.188 

1991 

326 

+9.814 

1991 

227 

+8J23 

1991 

329 

+10.004 

1991 

232 

+8.665 

1991 

336 

+10.449 
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Table  5.  Timing  Pulse  Differences  (con.) 


YEAR 

DAY 

DIFFERENCE 

YEAR 

DAY 

DIFFERENCE 

1991 

340 

+10.701 

1992 

073 

+1.704 

1991 

343 

+10.892 

1992 

076 

+1.875 

1991 

347 

+11.140 

1992 

080 

+2.095 

1991 

354 

+11.573 

1992 

083 

+2.244 

1992 

003 

+12.091 

1992 

087 

+2,441 

1992 

006 

+13.484 

1992 

092 

+2.692 

1992 

023 

+13.424 

1992 

097 

+2.940 

1992 

030 

+13.809 

1992 

106 

+3J95 

1992 

034 

+14.073 

1992 

113 

+3.754 

1992 

030 

+16.019 

1992 

120 

+4.112 

1992 

041 

+15.213 

1992 

125 

+4J90 

1992 

045 

+14.970 

1992 

153 

+6.179 

1992 

050 

+15.293 

1992 

157 

+6.440 

1992 

055 

+15,610 

1992 

160 

+6.650 

1992 

057 

RESET 

1992 

167 

+7.114 

1992 

059 

+0.853 

1992 

174 

+7.620 

1992 

062 

+1.040 

1992 

181 

+6.173 

1992 

066 

+1.602 

1992 

188 

+6.985 

1992 

069 

+1.469 

1992 

195 

+7.470 
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Table  5.  Timing  Pulse  Differences  (con.) 


YEAR 

DAY 

DIFFERENCE 

YEAR 

DAY 

DIFFERENCE 

1992 

199 

+7.567 

1993 

042 

+13.945 

1992 

209 

+7.760 

1993 

049 

+14J81 

1992 

213 

+7.615 

1993 

091 

+16.740 

1992 

220 

+7.913 

1993 

098 

+12.586 

1992 

232 

+5.688 

1993 

105 

+17.446 

1992 

241 

+6.279 

1993 

112 

+22.776 

1992 

244 

+6.472 

1993 

118 

+18.427 

1992 

250 

DOWN 

1993 

139 

+18.484 
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Ihiiveml  Time  (UT).  Identifying  these  time  sub-intervals  of  linear  POGS  clock-drift  requires 
careful  analysis  of  the  time  diffoence  data  in  Figure  4  and  the  comparison  of  that  data  with  the 
magnetic  total  intensity  data.  The  magnetic  field  data  exhibit  small  anomalies  when  the  POGS 
clodt  makes  a  discontinuous  break  from  its  linear  drift  trend.  In  this  way  the  start-time  and 
end-time  oi  each  linear  clock-drift  time  sub-interval  can  be  determined  to  within  a  few 
milliseconds.  However,  ionospheric  and  magneto^heric  magnetic  disturbances  can  occasionally 
interfere  with  and  reduce  the  accuracy  of  this  process. 

In  addition  to  the  clock-drift  corrections,  there  is  one  rather  small  time  correction  of  0.51 
milliseconds  applied  to  the  data.  This  correction  arises  from  a  timing  delay  on  the  POGS 
clock-timer  board.  When  all  of  the  timing  corrections  are  applied  to  the  data,  the  overall  timing 
error  associated  with  the  data  is  expected  to  be  on  the  order  of  10  to  20  milliseconds  for  most  of 
the  data.  The  timing  errors  of  the  ephemeris  data  are  of  the  same  magnitude.  Both  data  sets, 
having  been  synced  with  Naval  Observatory  time  to  the  same  accuracy,  are  then  said  to  be 
synced  with  each  other  and  are  then  merged  by  time.  Since  the  ephemeris  data  is  sent  to 
NAVOCEANO  at  1-minute  time  intervals,  it  must  be  interpolated  in  order  to  attach  a  position  to 
each  lO-second  magnetic  field  measurement.  Cubic  spline  interpolation  is  used  on  the 
Earth-Centered,  Earth-Fixed  rectangular  coordinates  of  the  DMA-supplied  ephemeris. 

6.  DATA  EDITING 

The  observed  magnetic  field  data  are  not  conq)letely  clean.  During  the  course  of  a 
single  orbit,  the  satellite  encounters  numerous  radiation  hits,  some  of  which  may  alter  conq)uter 
drip  manory  bits  through  some  ionization  process.  The  altered  data  appear  as  a  set  of  spikes  in 
the  total  intensity  as  indicated  in  Figure  S.  These  spikes  are  flagged  using  an  interactive 
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Figure  5.  FOGS  Uncalibrated/Unedited  Ibtal  Intoisity  Residual  Profile  Computed  with  Respect  to  the  WMM-90  Modd: 
Covers  13379.1  Seconds  (2  Orhits);  the  Large  Oscillations  are  the  Phase  Differences  Betwem  VOGS  Data 
Uncorrected  for  Clock  Drift  and  Ae  WMM-90  Model;  the  Spikes  are  Due  to  Radiation  Encountered  in  Space;  the 
Interval  Betwem  Small  Tick  Marks  on  the  ''i^ical  Axis  Corresponds  to  100  nT 


grqihics  editor.  An  edit  flag  Mt  to  "0"  for  good  data  and  set  to  T  for  bad  data  is  attached  to 
each  observation.  A  siinilar  flag  is  used  to  indicate  when  the  tracking  beacons  are  turned  on  and 
off.  A  data  recod  is  considered  to  be  a  valid  magnetic  field  measurement  when  both  flags  are 
zero.  * 


7.  PRELIMINARY  MODELING  OF  POGS  DATA 

The  POGS  satellite  oibit  is  a  polar  (me  and  has  an  average  altitude  of  t^proximately  688  km. 
The  most  significant  ionoifdmric  current  systems  are  in  the  altitude  range  of  100  km  to  200  km. 
These  current  systems  include  die  auroral  electrojets  in  the  north  and  soudi  polar  regions  and  the 
eqoatrnial  dectrojet  vdiich  straddles  the  ^cmaagnedc  equator.  Althou^  the  currents  in  these 
regions  are  generally  well  defined  and  Icmalind  excqit  for  periods  of  magnetic  storms,  the 
fidds  duU  diey  generate  fiequendy  affect  the  ambient  magnetic  field  at  mid-latitudes. 
Additionally,  diere  are  a  variety  of  magnetot^heric  currents  located  several  earth  radii  above  the 
POGS  satellite  cvbit  dut  also  generate  fidds  diat  can  be  sensed  by  the  POGS  satellite's 
Experiment  magnetometer.  The  most  significant  magnetoqiheric  source  of  magnetic  fields  is 
die  Ring  cunrent,  which  surrounds  die  eardi  in  die  equaUnial  plane.  Finally,  there  are 
field-aligned  current  systems  which  coiqile  the  magneto^here  to  the  iono^here  in  the  polar 
r^ions.  Since  POGS  is  a  polar  (nbiting  satellite,  it  will  pass  direcdy  through  the  field-aligned 
currents,  and  during  certain  periods  of  hdghtened  solar  activity,  its  magnetometer  may  measure 
magnedc  fidds  generated  by  diese  fidd-aligned  currents  on  the  order  of  a  thousand  nanoTedas 
or  more. 

It  is  our  intended  purpose  to  globally  modd  only  the  main  (core-goierated)  naagnedc  field  of 
die  eardi  using  POGS  data.  This  means  that  the  remnant  and  induced  magnetic  fields  generated 
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by  the  eutii't  crust,  as  well  as  those  magnetic  fields  generated  by  exoatmo^heric  electric 
currents  that  exist  between  the  iono^here  and  the  magneto^here,  must  be  separated  fit>m  those 
generated  within  the  earth's  core.  This  separation  is  accomplished  in  several  diffoent  ways. 
The  hei^t  of  the  satellite  above  the  eardi's  surface  is  an  effective  filter  of  the  crustal  fields.  To 
assist  in  the  effcnt  of  separating  the  exoatmopheric  fields,  the  National  Geophysical  Data  Cotter 
(NGDQ  in  Boulder,  CO,  computes  the  planetary  K  (K^  index  which  mnniitors  the  magnetic 
effect  of  solar  activity  at  3-hour  intervals.  This  index  ranges  in  value  from  1  to  10,  with  1 
coneqxmding  to  solar  quiet  periods  and  10  correponding  to  solar  disturbed  periods. 

Our  first  task  is  to  select  POGS  data  on  the  basis  of  the  K,  index.  Goierally,  data  wiOt  an 
index  equal  to  or  less  than  2*,  >i(diich  correpcmds  to  the  most  magnetically  quiet  times,  are 
selected  for  modeling.  Then,  from  this  magnetically  quiet  set  of  data,  we  select  data  on  the  basis 
of  local  time.  The  mu^etically  quietest  part  of  the  day  is  betweoi  9  p.m.  and  6  a.m.  We  also 
sdect  data  m  die  basis  of  season.  Goierally,  the  magnetically  quietest  part  of  the  year  is  during 
the  three  winter  iiKmths  for  the  north  and  south  hemipheres.  Additionally,  a  priori  global 
magnetic  field  models  do  exist  vdiich  can  be  used  as  refooices.  The  residual  fields  from  these  a 
priOTi  models  can  be  used  to  visually  edit  and  remove  from  the  data,  via  interactive  computo 
gnqdiic  techniques,  the  effect  of  the  exoatmopheric  currents  on  the  satellite  measurements. 
External  fidds  due  to  magnetophoic  and  ionopheric  effects  can  also  be  modeled  out  in  a 
plwmniienological  manner  using  the  Disturbance  Storm  Time  (Dst)  index,  wdiich  is  another 
measure  of  magnetic  activity,  usually  attributed  to  Ring  current  effects,  but  more  recendy  also 
attributed  in  large  part  to  field-aligned  current  effects.  An  example  of  the  Dst  index  generated  by 
M.  Sugiura  fw  the  Magsat  mission  (1979  to  1980)  is  illustrated  in  Figure  6.  For  solar  quiet 
times  in  dw  dawn-dusk  orbit  of  Magsat,  this  index  was  typically  less  than  50  nT  and  had  a 
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Figure  6.  The  Dst  Index  During  One  Magsat  Orbit  (Units:  nT) 


distiiictive  latitudinal  variation  as  shown.  A  revised  and  in^roved  method  of  generating  the  Dst 
index,  \riiich  we  call  the  Modified  Dst  index,  is  being  genmted  for  the  POGS  data  by  Wallace 
Campbell  (in  Press)  of  the  U.S.  Geological  Siwey  in  conjunction  with  NGDC.  During  quiet 
times,  the  magnitude  of  this  index  for  POGS  data  will  be  substantially  greater  than  50  nT  since 
POGS  will  see  all  local  time,  not  just  dawn  and  dusk  as  with  Magsat.  To  minimize  this  source 
of  noise,  we  restrict  our  data  selection  to  the  night  side  of  the  earth. 

Finally,  we  can,  as  part  of  the  core-field  spherical-harmonic  noodeling  process,  include  terms 
in  the  model,  which  correspond  to  external  field  effects  and  thereby  model  out  the  longer 
wavelength  magneto^heric  contributions  to  the  measured  field.  The  only  exoatmospheric 
mr^etic  fields  for  which  we  caruiot  a^rcount  for,  even  in  a  phenomenological  maimer,  are  the 
mid-latitude  ionospheric  curren*  system  effects  even  though  much  is  known  about  them 
(Campbell  [1989]).  Generally,  these  are  weak  current  systems,  particularly  during  magnetically 
quiet  times  for  which  the  index  is  2^  or  less,  and  for  the  most  part  are  ignored. 

The  model  presented  in  this  report  is  preliminary  in  nature.  This  means  that  the  POGS  data 
selected  for  this  model,  though  edited  for  radiation-induced  iqiikes  and  POGS  beacon 
transmissions,  were  not  edited  for  exoatmosphmc  current  effects.  The  data  were  selected  for  a 
index  of  2  or  less.  The  data  were  also  selected  to  reside  within  local  times  between  9  p.m. 
and  3  a.m.  in  order  to  minimize  the  effect  of  magnetic  fields  generated  by  solar 
heating/ionization  of  the  upper  atmosphere  and  ionosphere.  The  data  selected  covered  the 
mmths  of  January  and  February  1991.  Consequently,  the  winter  season  was  in  progress  in  the 
nofthem  hemi^here,  vriiile  the  summer  season  was  in  progress  in  the  southern  hemi^here.  No 
seastmal  ccwrections  were  performed.  Furthermore,  no  corrections  for  Dst  woe  poformed  since 
the  Dst  indices  were  not  yet  available. 
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7.1  The  POGSPrcilniBjiry  Model 

A  8|dierical-hannoiuc  model  through  internal  degree  12  and  extoiial  degree  2  was  generated 
using  40,491  total  intensity  POGS  data  values,  with  a  mean  Epoch  t  corresponding  to  day  34  of 
1991.  The  potential  field  model  is  of  the  following  mathematical  form: 

V(r,e,9,x)  =  a  2  £  (^J  co8m<p  +  h«(x)  8inm<|>]Pr(e)  + 

(8) 

a  X  2  (■!■)  [<!•-('')  cosnMP  +  8inm(|>]P?(e) 

■  a  I  BaO 

where  die  parameter  "a"  is  the  mean  radius  of  the  earth  (6371.2  km),  the  parametos  g^  and  h^ 
are  the  internal  Gauss  coefficients,  and  the  parameters  and  s^  are  the  external  Gauss 
coefficients.  The  modeling  objective  was  to  determine  these  Gauss  coefficients  at  the  mean 
q)och  X  of  die  data  set  Subsequendy,  the  internal  Gauss  coefficients  were  temporally  adjusted 
to  the  1990.0  Epodi  using  die  Secular  Variati(m  portion  of  die  1990.0  Epoch  International 
Geomagnetic  Reference  Field  (IGRF-90)  (lAGA  WG  V«8  [1992])  in  order  to  compare  the 
POGS  inetiminary  model  diiecdy  to  other  models  of  thatqioch. 

The  magnetic  field  B  is  the  negative  gradient  of  this  potential  field  V.  That  is: 

B(r,  e,  q>,  X)  *  -  V  V(r,  0, 9,  x)  (9) 

The  itdemal  Gauss  codficients  are  assumed  to  be  linearly  dependent  on  time  and  are 
considered  valid  for  tune  intervals  on  the  ruder  of  S  years  or  less.  These  coefficients  are 
generally  referenced  to  abase  qiochof  time T^^,  which  in  units  of  years  is  divisible  by  5. 
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Thus,  tfiey  take  the  following  mathematical  forms; 


gMi(^)  ®  gaalT Epoch)  "*■  £•<■  (X  ~  TEpoch)  (1®*) 

haaiCT)  s=  haaCr epoch)  +  h.*  (t  -  T^h)  (10b) 

The  coefficients  and  1^  have  units  of  nanoTeslas  as  do  the  coefficients  and  s^  .  while 

•  • 

and  h^  are  the  secular  rates  of  change  of  their  emresponding  internal  Gauss  coefficients. 

The  extanal  Gauss  coefficients  are  assumed  to  dq>end  on  the  time-dependent  Disturbance 
Suwm  Time  index  Dst(t),  when  it  becomes  available,  as  follows: 


-  Om  -f 

(Ha) 

s,.(t)  *  y,«  +  6«iD*(t) 

(lib) 

In  order  to  avoid  generatif^  ^Mirious  anomalies  into  the  model  via  the  Backus  effect  (Stem 
et  al.  [1980]),  vector  data  from  the  IGRF-90  model  were  foliated  at  l‘  intervals  throughout  a 
region  straddling  ±  2(f  about  the  geomagnetic  equator.  The  data  were  unwei^ted.  The 
resulttng  model  coefficients  are  given  in  Table  6,  where  they  are  compared  with  the  1990.0 
Epoch  U.'VUK  World  Magnetic  Model  (WMM-90)  coefficients,  which  were  based  on 
1979/19^  vint^  Magsat  data  and  on  a  Project  MAGNET  data  set  which  spanned  th«  decade 
of  the  1980's  and  which  has  a  mean  age  of  S  years  with  respect  to  the  1990  modding  qxxh 
(Quhm  et  al.  [1991]).  The  POGS  model  coefficients  listed  in  Table  6  were  temporally  adjusted 
to  die  1990  epoch  from  the  mean  modeling  epoch  of  the  POGS  data,  unng  the  IGRF-90  secular 
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Table  6.  Preliminary  POGS  Spherical  Harmonic  Model  Gauss  CoefOcients 
(IntemaO  Compared  to  WMM-90  at  1990.0;  Units:  nT 


a 

m 

EQGS 

WMM-90 

8m 

1 

0 

-29S00.8 

0.0 

-29780.5 

0.0 

1 

1 

-1852.8 

5413.2 

-1851.7 

5407.2 

2 

0 

-2113.5 

0.0 

-2134J 

0.0 

2 

1 

3040.5 

-2283J 

3062.2 

-22783 

2 

2 

1493.7 

-384.6 

1691.7 

-3843 

3 

0 

1305.1 

0.0 

1312.9 

0.0 

3 

1 

-2241.5 

-284.7 

-2244.7 

-284.9 

3 

2 

1245.0 

294.7 

1244.8 

291.7 

3 

3 

807.4 

-353.5 

808.4 

-352.4 

4 

0 

937.4 

0.0 

933.5 

0.0 

4 

1 

783.9 

247.4 

784.9 

249.4 

4 

2 

325.1 

•236.9 

323.5 

-232.7 

4 

3 

•419.4 

87.1 

-421.7 

913 

4 

4 

138.7 

-294.7 

139.2 

-294.5 

S 

0 

•204.1 

0.0 

-2083 

0.0 

S 

1 

354.7 

39J 

352.2 

40.8 

5 

2 

U53 

154.0 

254.5 

148.7 

S 

3 

-113.0 

-151.5 

-110.8 

-154.6 

S 

4 

-144.5 

-46.8 

-142J 

-47.6 

5 

5 

-37.2 

98J 

-37.2 

97.4 

4 

0 

58.7 

0.0 

59.0 

0.0 

4 

1 

40.9 

-9.0 

43.7 

-14.7 

4 

2 

41.9 

79.8 

40.0 

82.2 

4 

3 

-179.4 

44.8 

-1813 

70.0 

4 

4 

3J 

-57.7 

0.4 

-56.2 

4 

S 

17.5 

-1.4 

15.4 

-1.4 

4 

4 

-97.4 

23.0 

-94.0 

24.4 

7 

0 

75.9 

0.0 

74.1 

0.0 

7 

1 

-41.7 

-74.2 

-42.1 

-78.6 

7 

2 

2.5 

-24.0 

13 

-24.7 

7 

3 

27.9 

2.9 

30.2 

0.1 

7 

4 

OJ 

22.2 

4.7 

19.9 

7 

5 

4J 

18.8 

7.9 

17.9 

7 

4 

10.2 

-20.1 

10.1 

-21.5 

7 

7 

2.0 

-7.2 

1.9 

-4.8 
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Table  6.  Preliminary  POGS  Model  Comparison  (Con.) 


FOGS  WMM-90 


n 

m 

8mi 

g- 

8 

0 

25.8 

0.0 

22.9 

0.0 

8 

1 

2.2 

8.9 

2J 

9.7 

8 

2 

-3.5 

-19.7 

-1.2 

-19J 

8 

3 

-9.0 

6.9 

-11.7 

6.6 

8 

4 

-12.9 

-21.1 

-17.5 

-20.1 

8 

5 

5.5 

13J 

2.2 

13.4 

8 

6 

5.8 

10.1 

5.7 

9.8 

8 

7 

3.2 

-19J 

3.0 

-19.0 

8 

8 

-6.6 

-10.1 

-7.0 

-9.1 

9 

0 

5.6 

0.0 

3.6 

0.0 

9 

1 

7.4 

-19.4 

9.5 

-21.9 

9 

2 

3.2 

13.2 

-0.9 

14J 

9 

3 

-13.6 

11.2 

-10.7 

9.5 

9 

4 

7.8 

-6J 

10.7 

-6.7 

9 

5 

-4.8 

-7.4 

-3.2 

-6.4 

9 

6 

-1.8 

9.7 

-1.4 

9.1 

9 

7 

7.1 

7,7 

6J 

8.9 

9 

8 

OJ 

-8.1 

0.8 

-8.0 

9 

9 

-5.4 

2.2 

-5.5 

2.1 

10 

0 

-4.2 

0.0 

-3J 

0.0 

10 

1 

-3.0 

1.5 

-2.6 

2.6 

10 

2 

2.0 

0.6 

4.5 

1.2 

10 

3 

-4.9 

4.5 

-5.6 

2.6 

10 

4 

0.7 

4.1 

-3.6 

5.7 

10 

5 

5.8 

-4.6 

3.9 

-4.0 

10 

6 

1.8 

0.2 

3.2 

-0.4 

10 

7 

1.2 

-0.7 

1.7 

-1.7 

10 

8 

2.4 

3.2 

3.0 

3.8 

10 

9 

3.0 

-0.4 

3.7 

-0.8 

10 

10 

0.7 

-6J 

0.7 

-6.5 

11 

0 

2.6 

0.0 

1.3 

0.0 

11 

1 

•1.4 

-0.6 

-1.4 

0.0 

11 

2 

-IJ 

OJ 

-2,5 

1.0 

11 

3 

1.9 

-1.1 

3.2 

-1.6 

11 

4 

-0.5 

-2.4 

0.2 

-2.2 

11 

5 

-1.7 

0.2 

-1.1 

1.1 

11 

0 

-0.5 

-0.8 

OJ 

-0.7 

11 

7 

0.1 

-2.4 

-0.3 

-1.7 
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Table  6.  Preliminary  POGS  Model  Comparison  (Con.) 


n 

m 

K. 

WMM-90 

8» 

**B- 

11 

8 

OJ 

-1.8 

0.9 

-1.5 

11 

9 

-1.0 

-1.0 

-1.1 

-1.3 

11 

10 

2.5 

-1.2 

2.4 

-1.1 

11 

11 

3.0 

0.0 

3.0 

12 

0 

-1.8 

-IJ 

12 

1 

-1.7 

2.4 

0.1 

12 

2 

0.8 

1.2 

0.5 

12 

3 

1.0 

1.5 

0.7 

ij 

12 

4 

1.1 

-2.4 

0.4 

-1.5 

12 

5 

0.0 

0.1 

-0.2 

OJ 

12 

0 

-1.5 

0.8 

-1.1 

12 

7 

0.7 

-OJ 

0.9 

-1.1 

12 

8 

-0.0 

0.9 

-0.0 

1.2 

12 

9 

0.2 

0.1 

0.8 

12 

10 

03 

-0.9 

0.2 

-IJ 

12 

11 

0.5 

0.0 

0.4 

12 

12 

0.2 

0.0 

0.2 

0.0 
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variation  model  in  order  to  make  the  comparison  with  the  WMM-90  model  coefficients.  There 
are  significant  differences  between  some  comparable  lower  order  coefficients  of  the  two  models. 
This  is  partly  due  to  the  omission  of  ionospheric  and  magnetospheric  corrections  to  the  POGS 
data  and  partly  due  to  errors  in  the  WMM-90  model,  which  is  based  on  relatively  old  data. 

The  purpose  for  generating  this  preliminary  PCXjS  model  is  to  use  it  as  a  tool  for  further 
reduction  of  POGS  data.  In  particular,  the  magnetic  residuals  from  this  model  will  resolve  the 
field-aligned  current  effects  very  distinctly  and  thus  permit  the  removal  of  these  effects  in 
subsequent  models  based  on  this  data.  Thus,  the  process  of  generating  a  global  magnetic  field 
model  from  POGS  data  is  an  iterative  one.  The  model  in  Table  6  exhibited  the  following 
statistical  properties  with  respect  to  the  40,491  POGS  total  intensity  data  values  with  =  2^ 
collected  during  the  months  of  January  and  February  1991  that  were  used  to  generate  die  model; 

Mean:  -10.99  nT 

Standard  Deviation;  22.91  nT 

Root  Mean  Square;  25.41  nT 

These  statistics  include  the  contributions  of  both  the  internal  and  external  portions  of  the  POGS 
preliminary  model.  The  Gauss  coefficients  for  the  external  field  at  the  mean  epoch  of  the  POGS 
data  subset  are  listed  in  Table  7.  Since  the  preliminary  field  model  ignored  the  effects  of  the 
field-aligned  curroits  and  the  Dst  index  (i.e.,  =  0,  and  =0  for  all  n  and  m),  it  is  expected 

diat  subsequent,  more  detailed  models  will  reduce  these  statistics  substantially.  Nevertheless,  as 
they  stand,  the  statistics  are  good. 
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7.2  Analysis  of  POGS  Residuals 

The  residual  magnetic  field  is  composed  of  magnetic  fields  generated  by  all  sources  except 
the  earth's  outer  core.  There  is  a  tendency  to  expect  the  POGS  residuals  to  be  similar  to  Magsat 
residuals.  In  particular  it  is  often  commented  that  POGS  data  are  rather  noisy,  with  the 
inq>lication  that  there  may  be  something  wrong  with  the  calibration  of  the  magnetometer  or  with 
the  data  reduction  procedures.  Though  not  explicitly  stated,  the  implication  is  that  the 
con^>arison  is  made  with  respect  to  the  Magsat  data  set,  with  which  there  is  great  familiarity. 
There  is  actually  nothing  amiss  with  either  the  calibrations  or  the  data  reduction  procedures. 
What  some  modelers  seem  to  forget  is  that  the  Magsat  mission  was  sun-synchronous,  the 
satellite  traversing  in  a  dawn-dusk  orbit.  This  is  the  least  magnetically  disturbed  of  all  the 
possible  orbits  that  one  could  choose.  The  POGS  mission,  on  the  other  hand,  is  not 
sun-synchronous  and  therefore  sees  all  local  times.  Most  particularly,  it  sees  the  dayside  local 
times  between  6  a.m.  and  6  p.m.,  which  are  more  magnetically  active  due  to  day-side  solar 
heating  effects.  Secondly,  the  altitude  of  the  POGS  satellite  is  roughly  twice  as  high  as  that  of 
Magsat.  Consequently,  the  field-aligned  currents  through  which  POG^  traverses  are  broader 
latitudinally  than  those  through  which  Magsat  traversed  due  to  the  divergence  of  the  main 
magnetic  field  lines  with  altitude.  This  makes  the  job  of  removing  the  magnetic  effects 
geno^ted  by  field-aligned  currents  from  the  POGS  data  more  difficult  than  for  Magsat.  The 
problem  is  further  aggravated  by  the  fact  that  only  total  intensity  data  is  available  from  POGS, 
while  Magsat  had  vector  data.  Figure  7  illustrates  the  anomalous  magnetic  field  for  three  vector 
conqxments  observed  by  Magsat  as  it  traversed  through  the  freld-aligned  currents  feeding  the 
north  and  south  pole  auroral  ovals.  In  contrast,  a  comparable  POGS  orbit  is  exhibited  in 
Figure  8,  where  the  top  illustration  exhibits  total  intensity  data  collected  during  a  period  of 
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Vector  Magnetic  Field  Residuals  During  One  Magsat  Orbit  Computed 
with  Respect  to  the  GSFC 12/83  Model:  Clearly  Illustrates  the  Field-Aligned 
Current  Effect  in  the  Polar  Auroral  Zones;  the  Vertical  Axis  Range  is  ±  1000  nT 


modente  solar  activity  for  which  =  6.  The  magnitude  of  the  field-aligned  curroit  effect  is 
much  less  than  that  of  Figure  7  due  to  the  divergence  of  the  current  at  higher  altitudes.  For  the 
same  reason  diere  is  a  distinct  broadening  of  the  inner  and  outer  rings  of  the  ovals  in  die  two 
hemi^heres  and  a  genoal  spreading  of  the  oval  toward  the  equator.  The  second  graph  on 
Figure  8  illustrates  the  tenperatures  of  the  magneUHueter  sensor  electronics  box  (blue)  and  of 
the  magnetometer  sensor  (yellow).  The  tenperatures  are  characteristically  stable,  varying 
slowly  within  narrow  limits.  The  third  gnph  in  Figure  8  exhibits  the  satellite  attitude  relative  to 
the  sun-line.  There  is  no  attitude  information  when  the  satellite  is  on  the  night-side  of  the  earth. 

Figure  9  exhibits  data  collected  on  orbit  6113  during  a  highly  disturbed  period  for  vdiich 

==  8.  In  tlds  figure,  it  is  possible  to  see  a  broadening  of  the  auroral  ovals  in  both  hemispheres 
to  mid-latitudes  as  well  as  the  pronounced  effect  of  the  equatorial  electrojet  which  has  been 
enhanced  dramatically  by  magnetic  storm  activity.  The  total  intraisity  data  in  Figure  10  reflect 
modoate  solar  activity  corresponding  to  =  4^  to  5*.  Notice  that  at  the  equator  in  die  top  graph 
there  is  a  significant  difference  in  the  residuals  between  the  day-side  and  the  night-side  of  the 
earth  as  defined  by  the  attitude  data  in  the  third  graph  in  Figure  10.  On  the  day-side,  the 
equatorial  electrojet  contributes  an  additional  75  nT  to  100  nT  to  the  observed  field.  In  contrast. 
Figure  1 1  depicts  an  unusually  quiet  magnetic  field  record  collected  from  orbit  4106  for  which 
die  index  was  between  T  and  2'. 

A  low  index  does  not  always  indicate  benign  magnetic  field  effects  from  exoatmo^hoic 
sources.  Figure  12  for  orbit  4909  was  coll^ned  during  a  period  for  which  =  1.  At 
mid-latitudes,  particularly  in  the  southern  honisphere,  thore  is  a  voy  pronounced  anomaly 
v^ch  appears  to  be  caused  by  some  unexpectedly  strong  mid-latitude  current  system.  The 
residual  field  displayed  in  Figure  12  is  taken  with  respect  to  our  preliminary  FOGS  model.  The 
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Figure  8.  Total  Intensity  Residual  Profile  Computed  with  Respect  to  the  Preliminary 
POGS  Model  and  the  Temperature  and  Attitude  Profiles  for  Orbit  #5825: 


Figure  9.  Total  Intensity  Residual  Profile  Computed  with  Respect  to  the  Preliminary 

POGS  Model  and  the  Temperature  and  Attitude  Profiles  for  Orbit  #6113:  Kp  == 


Figure  10.  Total  Intensity  Residual  Profile  Computed  with  Respect  to  the  Preliminary 

POGS  Model  and  the  Temperature  and  Attitude  Profiles  for  Orbit  #6679:  =  4^  to  5* 
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Figure  11.  Total  Intensity  Residual  Profile  Computed  with  Respect  to  the  Preliminary 

POGS  Model  and  the  Temperature  and  Attitude  Profiles  for  Orbit  #4106:  K,,  =  1^  to  2' 


Figure  12.  Total  Intensity  Residual  Profile  Computed  with  Respect  to  the  Preliminary 

POGS  Model  and  the  Temperature  and  Attitude  Profiles  for  Orbit  #4909:  = 


anomalous  feature,  however,  is  independent  of  the  model  used.  The  magnetic  field  total 
intensity  residual  generated  from  the  IGRF-90  model  for  calendar  day  73,  1991,  exhibits  the 
same  feature.  Therefore,  this  anomaly  appears  to  be  in  the  data  and  seems  to  be  real,  but  is 
probably  of  external  origin.  Many  of  the  POGS  satellite  profiles  have  this  mid-latitude  feature, 
though  not  always  as  pronounced  as  is  illustrated  in  Figure  12.  Furthermore,  this  feature  persists 
in  the  southern  hemisphere  regardless  of  the  season.  There  appears  to  be  no  comparable  anomaly 
in  the  northern  hemisphere  regardless  of  the  season.  Until  additional  ionospheric  and 
magnetospheric  corrections  are  made  to  the  data,  our  initial  interpretation  of  this  anomaly  will 
remain  tentative. 

Finally,  a  direct  comparison  between  residuals  computed  from  the  POGS  prelinainary  model 
and  the  degree  10  IGRF-90  model  are  displayed  in  Figures  13  through  16.  Figure  13  illustrates 
the  POGS  residual  total  intensity  with  respect  to  the  POGS  preliminary  model  for  orbit  4044. 
Figure  14  is  the  IGRF-90  total  intensity  residual  for  the  same  orbit.  For  this  orbit  -  1.  The 
southern  hemisphere  mid-latitude  anonoaly  is  present  to  some  extent  in  both  figures.  Note  that 
on  1991  day  14  when  these  data  were  collected,  it  was  mid-winter  in  the  northern  hemisphere 
and  mid-summer  in  the  southern  hemisph^e.  As  a  second  contrasting  comparison,  an  orbit 
from  1991  day  173  was  selected,  which  corresponds  to  early  summer  in  the  northern  hemi^here 
early  winter  in  the  southern  hemisphere.  Figure  IS  illustrates  the  total  intensity  residuals  for  the 
POGS  preliminary  model,  while  Figure  16  illustrates  the  IGRF-90  model  total  intensity 
residuals  for  the  same  orbit.  The  mid-latitude  anomaly  is  not  as  evident  in  these  figures.  The 
index  for  6351  ranged  between  3*  and  4',  corresponding  to  moderately  active  external 
geomagnetic  activity  as  is  indicated  by  the  substantial  auroral  zone  field-aligned  current  activity. 
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Figure  13.  Total  Intensity  Residual  Profile  Computed  with  Respect  to  the  Preliminary 

POGS  Model  and  the  Temperature  and  Attitude  Profiles  for  Orbit  #4044:  Kp  - 


Figure  14.  Total  Intensity  Residual  Profile  Computed  with  Respect  to  the  IGRF-90 
Model  and  the  Temperature  and  Attitude  Profiles  for  Orbit  #4044: 


Figure  15.  Total  Intensity  Residual  Profile  Computed  with  Respect  to  the  Preliminary 

POGS  Model  and  the  Temperature  and  Attitude  Profiles  for  Orbit  #6351:  Kp  =  3^  to  4* 
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Figure  16.  Total  Intensity  Residual  Profile  Computed  with  Respect  to  the  IGRF-90 

Model  and  the  Temperature  and  Attitude  Profiles  for  Orbit  #6351:  -  3^  to  4' 


«.  FUTURE  PLANS  FOR  THE  1995  EPOCH  AND  BEYOND 


POGS  data  are  intended  to  support  the  199S  Epoch  World  Magnetic  Model.  In  the  coming 
months,  these  data  will  be  meticulously  edited  for  exoatmospheric  current  effects,  and  special 
Dst  corrections  will  be  implied.  Rather  than  using  data  with  indices  less  that  2* ,  as  was  the 
case  with  the  preliminary  model  and  other  models  based  on  Magsat  data,  only  data  with 
indices  less  than  T  will  be  used.  This  is  because  of  the  large  amount  of  POGS  data  that  are 
available.  Also,  only  data  from  the  nightside  and  only  data  from  the  winter  hemispheres  will  be 
used  since  seasonal  effects  do  seem  to  be  significant. 

In  die  preliminary  model  presented  in  this  report,  the  Baclcus  effect  was  controlled  by 
throwing  in  some  vector  data  from  a  predictive  model  along  the  geomagnetic  equator.  In  the 
1995  Epoch  model,  real  vector  data  collected  from  the  Project  MAGNET  aircraft  will  be  used. 
Higjh-level  Project  MAGNET  aeromagnetic  surveys  are  currently  being  executed  along  the 
geomagnetic  equator  in  the  Atlantic,  Pacific,  and  Indian  oceans.  Survey  tracks  that  are  either 
planned  for  the  next  6  months  or  have  already  been  collected  in  the  Atlantic  and  Pacific  oceans 
are  illustrated  in  Figures  17  and  18.  These  data  and  some  additional  data  collected  in  the  Indian 
ocean  will  be  sufficioit  to  control  the  Backus  effect. 

Because  of  the  longevity  of  the  POGS  satellite  coverage  and  because  of  the  spatial  density  of 
POGS  survey  coverage,  the  possibility  arises  that  a  sequence  of  models  can  be  made  to  perhaps 
degree  20  or  so,  at  6-month  intervals  covering  the  period  between  1991  throu^  1993.  This 
sequence  of  main  field  models  could  be  used  to  determine  for  the  geomagnetic  Secular  Variation 
model  for  that  time  period  to  degree  12  or  15.  This  is  in  contrast  to  those  of  only  degree  8  that 
are  currently  goioated  using  data  from  a  poorly  distributed  set  of  geomagnetic  observatories. 
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to  Counter  the  Backus  Effect  for  the  1995  Epoch  World  Magnetic  Model 


Finally,  with  respect  to  efforts  to  secure  data  for  future  modeling  epochs,  it  is  now  reasonably 
certain  that  there  will  be  a  series  of  boom-mounted,  scalar-quality,  vector-magnetic  experiments 
on  the  Block  5  DMSP  satellites,  which  will  span  the  time  fi^e  from  1995  through  2005.  Two 
of  these  satellites  will  be  maintained  in  orthogonal,  sun-synchronous  polar  orbits  at  all  times. 
Follow-on  Block  6  satellites  with  full  vector  capability  are  also  being  considered  to  cover  the 
2005  through  2015  time  frame.  Second,  we  note  that  the  Johns  Hopkins  Applied  Physics 
Laboratory's  Upper  Atmosphere  Research  Satellite  (UARS)  was  launched  in  September  1991, 
while  the  Swedish  satellite,  Freja,  was  launched  in  September  1992.  Attempts  are  being  made  to 
secure  vector  data  from  these  satellites,  although  the  data  reduction  may  be  somewhat  lengthy 
since  neither  data  set  is  of  Magsat  quality.  A  Danish  satellite  project  called  Oersted,  which  is 
intended  to  be  of  Magsat  quality,  had  been  planned  for  launch  in  the  latter  part  of  the  1990's. 
This  project  received  funding  from  the  Danish  government  in  1993  and  is  now  proceeding 
beyond  the  planning  stages  to  an  eventual  launch  in  cooperation  with  NASA  in  the  1996  time 
frame. 

The  current  prospects  for  securing  geomagnetic  survey  data  to  support  World  Magnetic 
Modeling  for  the  next  20  years  or  so  are  summarized  in  Figure  19.  It  seems  from  this  summary 
that  there  will  be,  as  this  century  closes  and  the  new  century  begins,  a  continuous  data  flow 
which  will  guarantee  the  availability  of  high-quality  Main  field  and  Secular  Variation  models  for 
the  foreseeable  future.  These  data  will  provide  long-term,  continuous  detailed  mapping  of  the 
earth's  magnetic  field  which  in  turn  will  provide  for  the  first  time  a  means  to  extend  our 
knowledge  of  the  earth's  deep  interior  through  close  monitoring  of  its  magnetic  secular  behavior. 
Thereby,  the  future  seems  bright  for  the  improved  predictive  capability  of  our  models  and  for 
the  wide  range  of  applications  they  support. 
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